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On the complex disc— corona interactions in the soft 
spectral states of Soft X-ray Transients 
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ABSTRACT 

Accreting black holes show a complex and diverse behaviour in their soft spectral 
states. Although these spectra are dominated by a soft, thermal component which 
almost certainly arises from an accretion disc, there is also a hard X-ray tail indicating 
that some fraction of the accretion power is instead dissipated in hot, optically thin 
coronal material. During such states, best observed in the early outburst of Soft X-ray 
Transients, the ratio of power dissipated in the hot corona to that in the disc can vary 
from ~ (pure disc accretion) to ~ 1 (equal power in each). Here we present results 
of spectral analyses of a number of sources, demonstrating the presence of complex 
features in their energy spectra. Our main findings are: (1) the soft components are not 
properly described by a thermal emission from accretion discs: they are appreciably 
broader than can be described by disc blackbody models even including relativistic 
effects, and (2) the spectral features near 5-9 keV commonly seen in such spectra can 
be well described by reprocessing of hard X-rays by optically thick, highly ionized, 
relativistically moving plasma. 
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1 INTRODUCTION 

Black Hole Soft X-ray Transient (SXT) sources are a sub- 
class of Low Mass X-ray Binaries occasionally undergo- 
ing dramatic outbursts. From a dim quiescent state they 
brighten (in the opt/UV/X/7-ray bands) by several orders 
of magnitude in the course of a few days, and then de- 
cline, with roughly exponential dependence of X-ray flux on 
time (see Tanaka & Lewin 1995; Tanaka & Shibazaki 1996; 
Chen, Shrader & Livio 1997 for reviews). The outbursts are 
thought to be caused by a sudden increase of mass accretion 
rate onto the central object as a result of the disc instability 
due to hydrogen ionization. The basic mechanism is thus the 
same as in dwarf novae (e.g. Osaki 1996), but with impor- 
tant differences due to larger and more massive accretion 
discs and the intense X-ray irradiation in SXT (Cannizzo 
1993; King & Ritter 1998). 

Black Hole SXT are perfect laboratories to study the ac- 
cretion process. Black holes are simpler accretors than neu- 
tron stars or white dwarfs thanks to the lack of global mag- 
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netic field and a hard surface. During an outburst and follow- 
ing decline the mass accretion rate changes by several orders 
of magnitude, giving data on how the accretion processes 
change as a function of rh. Results from SXT studies can 
then be applied to other systems with accretion onto black 
holes, e.g. active galactic nuclei (AGN). If similar outbursts 
take place in AGN, they cannot be observed in progress 
due to much longer (factor of 10 5 -10 7 ) time scales (Bur- 
den, King, & Szuszkiewicz 1998; Siemiginowska, Czerny, & 
Kostyunin 1996). 

Previous investigations of both black hole and neutron 
star SXT led to identification of a number of distinct spec- 
tral/temporal states (see e.g. Tanaka & Lewin 1995; van der 
Klis 1995 for reviews): (1) quiescence state, when the sources 
are very dim (Lx ~ 10 33 erg s _1 ), with (poorly constrained) 
power law spectra, (2) low/hard state, (LS) characterized by 
a hard, roughly power law spectrum (photon index r ~ 1.7) 
and a strong X-ray variability (r.m.s. 20-30 per cent), with 
power spectral density (PSD) peaking at ~ 0.1 Hz; (3) inter- 
mediate state (IS) , when the power law component is some- 
what softer (r ~ 2) and a weak soft, thermal component ap- 
pears in the spectrum; (4) high/soft (HS) state with X-ray 
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spectra dominated by a thermal component of temperature 
~ 1 keV and PSD of rather low amplitude (a few per cent) 
and a power law shape; (5) very high state (VHS) when 
both the soft, thermal (kT ~ 1 keV) and the hard power 
law components (T = 2 — 3) are present in the spectrum, 
and similarly the PSD is a mixture of PSDs characteristic 
to LS and HS. Additionally, quasi-periodic oscillations of 
various frequency (1-10 Hz) and strength appear in PSD in 
this state. For the purpose of this work the IS, HS and VHS 
are referred to as soft states, since it is in these states that 
a strong soft, thermal component is present in the spectra. 
Both the directly observed fluxes and temporal sequence of 
occurrence of these states during SXT's decline leave little 
doubt that the accretion rate is lowest in quiescence and 
increasing through LS, IS, HS to VHS. 

A general scenario for the evolution of accretion flow 
in SXT as a function of accretion rate was presented by 
Esin, McClintock & Narayan (1997). They suggested that 
the most characteristic feature of the evolution during the 
decline - the high to low state transition - involves a ma- 
jor change of the flow geometry. An optically thick, ther- 
mally emitting accretion disc dominant in soft states is re- 
placed by a hot, optically thin(ish) flow able to produce the 
strong, hard X-ray emission seen in the low state by thermal 
Comptonization (e.g. Gierlinski et al. 1997). The transition 
presumably involves evaporation of the optically thick disc 
(e.g. Meyer & Meyer-Hofmeister 1994; Rozariska & Czerny 
2000), whose inner radius increases from the last stable or- 
bit at 6 R g in the high state to ~ 10 4 7? g in the low state 
(R B = GM/c 2 is the gravitational radius). 

One way to track the optically thick accretion disc is 
through reflection of the hard X-ray spectrum by the accre- 
tion disc. The reflected component consists of the Compton 
reflected continuum with spectral features due to iron K 
shell absorption and fluorescence imprinted on it (Lightman 
& White 1988; George & Fabian 1991; Matt, Perola & Piro 
1991). When the reprocessing matter forms a disc rotating 
in a deep potential well, the spectral features are broadened 
and smeared by Doppler effects and gravitational redshift 
(Fabian et al. 1989; Ross, Fabian & Brandt 1996). The ob- 
served spectral features between 6-9 keV can be fit by such 
models and the derived amount of reflection and relativis- 
ts smearing is consistent with the disc inner edge moving 
inwards as the LS makes a transition to the IS/HS/VHS 
(Zycki, Done & Smith 1998; Gilfanov, Churazov & Revnivt- 
sev 2000). However, the amount of reflection and smearing 
seen in the LS generally precludes the disc being as far away 
as 10 4 R g as in the simplest versions of the truncated disc 
models (Esin et al., 1997), and instead requires that the 
optically thick material extends down to ~ 30 R s (Zycki, 
Done & Smith 1997; 1998; 1999a; Miller et al. 2001). While 
a (somewhat modified) truncated disc model can then fit the 
LS data, this is not a unique interpretation: the disc could 
extend down to the last stable orbit if there is complex ion- 
ization structure in the disc and/or relativistic outflow of 
the hard X-ray emitting plasma (see Done 2001 for review 
and the discussion in Di Salvo et al. 2001). 

At higher luminosities, in the IS/HS/VHS the spectrum 
is dominated by the soft component from the accretion disc 
(e.g. Ebisawa et al. 1994) as expected in the Esin et al. (1997) 
model, but this co-exists (especially in the VHS) with a hard 
X-ray tail. This tail is probably produced by non-thermal 



Comptonization (Coppi 1999; Gierlinski et al. 1999) as it 
can extend without a break to a few hundred keV (Grove 
et al. 1998). The ratio of luminosities in the hard and soft 
components can assume a broad range of values, from ~ 
where the spectrum is pure disc emission to ~ 1 with equal 
luminosities in the disc and X-ray hot phase (Ebisawa et al. 
1994; Homan et al. 2001; Rutledge et al. 1999; Sobczak et 
al. 2000) . As well as the continuum, iron spectral features at 
5-10 keV are again usually observed in the spectra (Ebisawa 
et al. 1994). They are interpreted as indicating reprocessing 
of hard X-ray radiation by an optically thick accretion disc, 
but detailed reflection models are not often fit to the data. 
A common spectral model is of a disc blackbody (Mitsuda 
ct al. 1984) and power law continuum with a broad iron line 
and smeared edge to phenomenologically model the reflected 
spectrum from the disc. While this model gives a rough idea 
of the ratio of soft to hard luminosity, it has some obvious 
drawbacks. Firstly, the disc spectrum presumably provides 
the seed photons for the Compton upscattering, so the hard 
spectrum cannot extend as an unbroken power law into the 
disc blackbody spectrum. A power law continuum overes- 
timates the contribution of the hard continuum at soft en- 
ergies, and so this suppresses and distorts the derived soft 
component. Also, a true reflected spectrum contains a con- 
tinuum as well as spectral features. Neglecting the contin- 
uum distorts the derived power law spectral index. The few 
papers which have modelled the reflected spectrum show 
unambiguously that the reprocessing medium in SXT in the 
IS/HS and VHS is highly ionized (Zycki et al. 1998; Miller 
et al. 2001; Wilson & Done 2001), unlike the LS where it is 
mainly neutral. Smearing of the reflected spectral features 
is observed in all spectral states (Zycki et al. 1997, 1998, 
1999a; Miller et al. 2001; Wilson & Done 2001). 

If the physics of the accretion flow really is controlled 
predominantly by m, then the SXT spectra should be the 
same as those in the persistent X-ray sources at similar ac- 
cretion rates. The most thoroughly studied soft state so far 
- in Cyg X-l (Gierlinski et al. 1999) - shows all the above 
listed features: non-thermal Comptonization as the origin of 
the hard X-rays, and the ionized and smeared reprocessed 
component. Additionally, it conclusively shows that the con- 
tinuum is not well modelled by a single disc blackbody spec- 
trum which is Comptonized by a non-thermal electron dis- 
tribution. Additional low temperature Comptonization of 
the disc spectrum is required by the data (Gierlinski et al. 
1999; Coppi 1999). 

In this paper we re-analyse some of the archival Ginga 
and RXTE data of SXTs in the soft states. We show that in 
most cases the spectra cannot be described by a physically 
motivated model of a simple (disc)blackbody emission, with 
a Comptonized power law tail which is reflected by the accre- 
tion disc. Similarly to the HS Cyg X-l spectrum described 
above, the soft component requires a broader spectral form 
than that of a disc blackbody (even including relativistic 
corrections). This could be due to additional thermal Comp- 
tonization in the disc, or to a separate blackbody component 
from X-ray heated hot spots on the disc. The poor fit of the 
disc blackbody model is contrary to most results derived 
from phenomenological models of the hard component and 
spectral features (power law plus broad line and/or smeared 
edge) which generally give adequate fits to the data, and 
on the few occasions where they did not (e.g. Ebisawa et 
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al. 1994), then the evident problems in the modeling made 
any inference about the soft component shape unconvinc- 
ing. The phenomenological models for the spectral features 
have more free parameters than the physically motivated 
reflected spectral models, and so can compensate for diffi- 
culties in the continuum model. We demonstrate that when 
the spectral features at 6-9 keV are modelled as X-ray re- 
processing, the reprocessor is strongly ionized, as expected 
for a disc at temperature > 0.5 keV. 



2 DATA AND MODELS 

We analyse Ginga Large Area Counter (LAC) data of 
GS 1124-68 (Nova Muscae 1991) and GS 2000+25. Gmga 
data are still some of the best to study spectra with strong 
soft components, thanks to its broad band energy coverage 
from ~ 1 — 20 keV. The data were re-extracted from LED AS 
public archive at Leicester University. The background was 
subtracted according to the 'universal' procedure based on 
blank sky observations (Hayashida et al. 1989). Whenever 
the source counts contaminated the Surplus above the Upper 
Discriminator (SUD) monitor, we first recovered the origi- 
nal SUD values using procedure described in Zycki et al. 
(1999a), and then applied the 'universal' background sub- 
traction. We assume 0.5 per cent systematic errors in the 
data. This is motivated by the uncertainty in LAC response 
from fits to the Crab spectrum (Turner et al. 1989), as well 
as the residual uncertainty in the determination of the LAC 
background (Hayashida et al. 1989) . 

We note that for data from instruments like Gtnga/LAC 
or RXTE/PCA the systematic error is often the dominant 
error up to about 10 keV. Its value then has a direct in- 
fluence on the xl values obtained. While this makes the 
absolute value of xi and the confidence contours on fit pa- 
rameters lose their usual statistical meaning, the \v values 
do preserve the ranking of models. Since some of our best 
fits have xl substantially less than 1, it suggests that the 
assumed systematic error may be somewhat overestimated. 

The RXTE data (PCA top layer only, detectors 0, 
1 and 2 only) were extracted from public archive at 
HEASARC/GSFC, and reduced with FTOOLS ver. 4.2 
software. The background was estimated employing the lat- 
est background models appropriate for bright sources (the 
'sky-VLE' model for Epoch 3 of RXTE observations). Sys- 
tematic errors of 1 per cent were assumed in each data chan- 
nel. 

The Ginga data we model here were previously anal- 
ysed by Ebisawa (1991) and Ebisawa et al. (1994). They 
used the phenomenological model of a disc blackbody (Mit- 
suda et al. 1984) plus a power law with smeared edge 
and broad gaussian line near 6-8 keV. Here we model the 
spectral features using the rel-repr model (Zycki et al. 
1999a). The model uses the angle-dependent Green's func- 
tions for the problem of Compton reflection as computed 
by Magdziarz & Zdziarski (1995; as implemented in XSPEC 
models pexrav/pexriv), to find the reflected continuum for 
a given primary continuum. The Fe K a line is computed as in 
Zycki & Czerny (1994), with photo- ionization computations 
as in Done et al. (1992), and is then added to the reflected 
continuum, so that the properties of the reflected continuum 
and the line are computed consistently for a given ioniza- 



tion state, inclination and elemental abundances. Relativis- 
tic and kinematic corrections are then applied by convolving 
the spectrum with the the appropriate Green's function for 
Schwarzschild metric, using the prescription of Fabian et 
al. (1989). The final model thus have three main parame- 
ters: (1) the ionization parameter, £ = 47rF/nH, where F is 
the irradiating flux in 5 eV - 19 keV band and tih is the 
hydrogen number density, (2) the overall amplitude (nor- 
malized to the underlying continuum), which is equivalent 
to the solid angle of the reprocessor as visible from the X- 
ray source normalized to 27r, fi/(27r) and (3) the inner disc 
radius, 7?i n determining the relativistic smearing of the spec- 
tral features (the irradiation emissivity is assumed oc r -3 , 
Zycki et al. 1999a). All radii are given in units of the gravi- 
tational radius (R s = GM/c 2 ). 

We note that pexriv and rel-repr are based on com- 
putations of irradiated discs with constant density. In more 
advanced computations with hydrostatic equilibrium, the 
thermal instability of X-ray irradiated plasma plays a role, 
and leads to a more complex vertical structure, with impor- 
tant observational consequences (Rozariska & Czerny 1996; 
Nayakshin, Kazanas & Kallman 2000). However, the ther- 
mal instability is probably unimportant for the soft state 
spectra in SXT's as the thermal equilibrium curve (the ion- 
ization parameter - temperature, S—T diagram) degenerates 
since its lower stable branch is at high disc temperature 
To ~ 0.5 — 1 keV. The soft spectrum dominates the spec- 
trum, so the Compton temperature is not much higher than 
To and the instability is suppressed. 

We use two analytical models approximating the 
inverse-Compton spectra: thComp, based on solution of 
the Kompaneets equation (Zdziarski, Johnson & Magdziarz 
1996) and COMPTT (Titarchuk 1994). We check whether our 
conclusions may be affected by the approximate character 
of the Comptonization models, using the COMpPS Comp- 
tonization model (Poutanen & Svensson 1996), which finds 
a numerical solution of the Comptonization problem for op- 
tically thin (r < 1) plasma explicitly considering succes- 
sive scattering orders. This code can be used with a variety 
of electron energy distributions, including a purely thermal 
distribution, a purely non-thermal one or a hybrid electron 
distribution, i.e. a Maxwellian distribution with a power law 
tail above certain Lorenz factor, 71 (Coppi 1999). Such a hy- 
brid electron distribution can provide a good description of 
X-ray/7-ray spectra of Cyg X-l in the soft state (Gierliriski 
et al. 1999) as Compton scattering of the disc blackbody by 
the thermal part of the electron distribution provides an ad- 
ditional, broad soft component while the X-ray power law 
is produced by Compton scattering of the disc photons by 
the non-thermal electron tail. 

The soft component will be modelled by the diskbb 
model (Mitsuda et al. 1984). We will also use a general rel- 
ativistic disc spectrum model GRAD, which integrates black 
body spectra over disc radius, taking into account: the colour 
temperature correction (constant with radius) at T C oi/T e g = 
1.7 (Shimura & Takahara 1995), fully relativistic radial de- 
pendence of disc temperature in Schwarzschild metric (Kato, 
Fukue & Mineshige 1998) and kinematic/relativistic effects 
on photons propagation, using the method of Fabian et al. 
(1989). A similar model for the Kerr metric, with the black 
hole spin a = 0.998 will be also used. Here we employ the 
radial effective temperature prescription from Kato et al. 
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(1998), the colour temperature correction constant with ra- 
dius, while the photon propagation is done by convolving 
the local spectra with the transfer function of Laor (1991), 
as implemented in XSPEC. We note that the colour tem- 
perature correction to a blackbody spectrum is appropriate 
when the escaping photons undergo saturated Comptoniza- 
tion. 

We use XSPEC ver. 10 for all spectral fits (Arnaud 
1996), with the above non-standard models implemented 
as local models. Photoelectric absorption is modelled using 
WABS (Morrison & McCammon 1983), and parameter errors 
are given as A% 2 = 2.7. 



3 INTERMEDIATE STATE DATA 

First, we analyse a number of datasets of various sources 
obtained when the sources were about to make a transition 
from high state to low state. We refer to them as Interme- 
diate State data. 



3.1 GS 2000+25 

We first analyse the Ginga 1-20 keV spectra of GS 2000+25, 

beginning with the Dec 8, 1988 dataset (day 224 after the 

peak on April 28, the last observation in the HS; see Ebi- 

sawa 1991, Tanaka & Lewin 1995, fig. 3.11). Repeating the 

analysis of Ebisawa (1991), we use the diskbb model and 

a power law tail with the smeared edge applied only to the 

latter component, wabs*(diskbb + smedge*(powerlaw)) 

(model A in Table |l|). This phenomenological description 

gives a rather poor fit, xt = 38/23, with best fit Ah = 
q+0.7 x 1Q 2i 

cm , i.e. lower than the likely interstellar value 
of ~ 6 x 10 21 cm" 2 . Adding a Gaussian line to the model 
does not improve the fit irrespectively of whether the line 
is narrow or broad, and its energy fixed at 6.4 keV or 
free. Replacing the disc blackbody model by the Comp- 
tonized blackbody (thComp) results in significant improve- 
ment, xt = 17.9/23 (with Ah fixed at the interstellar value; 
model B in Table [l]). The fit is even better when the hard 
power law with smeared edge is replaced by Comptoniza- 
tion model thComp (so that the full model contains two 
thComp components) with its corresponding reprocessed 
component, REL-REPR: the fit has xt — 10.6/23 (model D). 
As a test, we return to the diskbb model for the soft com- 
ponent, but we keep the thComp+rel-repr model for the 
hard component (model C in Table |l|) . We obtain a signifi- 
cantly worse fit, with xt = 30.0/23, even though Ah was free 
in this fit (x 2 increases to 99.7 when Ah is fixed at the inter- 
stellar value). Thus the best fit is obtained with the proper 
physical description of the Fe spectral features and a Comp- 
tonized blackbody model for the soft component (model D 
in Figure [jj . 

A very similar sequence of results is obtained for the 
earlier observations of GS 2000+25 (Nov 5th, Oct 18th, Sep 
8th, see Table ^), when the source luminosity was higher 
by up to a factor of ~ 5 than on Dec 8th. In all cases bet- 
ter fits are obtained with the soft component modelled as 
Comptonized blackbody than a disc blackbody, irrespective 
of the model used for the hard component (a power law with 
smeared edge or thComp with rel-repr, although the lat- 
ter is a more physical description of the source and gives a 
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Figure 1. Upper panel: The unfolded spectrum of GS 2000+25 on 
Dec 8 1988, shortly before its transition to the low state, showing 
a Comptonized soft component, Comptonized hard power law 
tail and its X-ray reflection from the disc (model D in Table 1). 
Lower panel: contributions to \ 2 from different models (Table 1). 
The best fit model contains Comptonized soft component and the 
physical description of X— ray reprocessing: thComp + rel-repr. 



better fit). Since the power law slope of the soft component 
is rather high (steep spectrum), the Comptonization tem- 
perature, kT e , is not well constrained by thComp. Thus the 
soft component can be well fit by either low temperature, 
high optical depth plasma or by high temperature, low opti- 
cal depth plasma. However, for the latter case when r c < 1, 
the approximations on which thComp is based become in- 
valid. Accurate Comptonized spectra are not a power law 
as the individual scattering orders can be seen, as can the 
initial seed photons. Using the COMpPS model (which gives 
accurate Comptonized spectra including the seed photons 
for r e < 1) for the soft component gives an upper limit on 
the temperature of ~ 25 keV, corresponding to r c « 0.55 
for spherical geometry for the Sep 8th spectrum. Thus these 
data favour a low temperature, optically thick(ish), unsat- 
urated (y = 4kT e /(m e c 2 ) r c ~ 0.1) Comptonization of the 
soft component. 
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Table 1. Results of model fitting of the Dec 8th spectrum (IS) of GS 2000+25 



parameter 


units 


A 


B 


C 


D 


Nn 


10 21 cm' 2 


0+0.6 


6.0(f) 


0+0.7 


6.0(f) 


T 
kT e 

r SO ft 


keV 
keV 


0.32 ±0.01 


0.18 ±0.1 

10(f) 
7.1 ± 0.3 


0.32 ±0.01 


0.18 ± 0.01 
10(f) 
7.3 ± 0.3 


-^edge 
7"cdgc 


keV 


8 - L -0.6 

0.52 ±0.28 


7.8 ±0.3 
32 ± 9 






r 




2.39 ±0.07 


1.93 ± 0.15 


2 54+ ' 07 


2 30+ - 11 


Q/2tt 

e 

Rin 


erg cm s — 1 

Rg 






0.44 ±0.32 
> 20 


U - JD -0.13 

(lt 3 ° g ) x 10 3 
> 8 


X v 




38/23 


17.9/23 


30.0/23 


10.6/23 



Model A: absorption (diskbb + smedge(powerlaw)) 
Model B: absorption(thComp + smedge(powerlaw)) 

Model C: absorption(diskbb + thComp + rel-repr) 
Model D: absorption(thComp + thComp + rel-repr) 



Table 2. Results of model fitting (x 2 /dof) of the Intermediate State spectra of GS 2000+25 



datasct 


model A 


model B 


model C 


model D 


model 


model E 


model F 


model G 


Sep 8th 


28.9/24 


11.0/23 


90.0/24^ 


14.2/24^ 


199/26 


15.6/25 / 


18.7/23 / 


66.7/24 


Oct 18th 


43.9/24 


14.5/23 


159/24 


17.2/23 


314/26 


21.0/24 


12.5/22 


80.7/24 


Nov 5th 


28.5/24 


14.5/23 


113/24 


14.4/23 


404/26 


19.2/24 


12.1/22 


60.6/24 


Dec 8th 


38.0/23^ 


17.8/22 


30.0/23^ 


10.6/23^ 


46.3/25^ 


15.1/25J" 


9.0/22-f 


32/23^ 



Model A: absorption(diskbb + smcdge(powerlaw)) 
Model B: absorption(thComp + smedge(powerlaw)) 
Model C: absorption(diskbb + thComp + rel-repr) 
Model D: absorption (thComp + thComp + rel-repr) 
Model 0: absorption(thComp + thComp), i.e. model D without the reprocessed component 
Model E: absorption(compPS + rel-repr), hybrid Comptonization 
Model F: absorption(diskbb + bbody + thComp + rel-repr) 

Model G: absorption(grad + thComp + rel-repr) 
N - best fit TVh significantly lower than the interstellar value 
' - Nn fixed at the interstellar value 



The COMpPS model can also be used with a hybrid 
electron distribution (both thermal and non-thermal elec- 
trons). This hybrid Comptonization model provides a good 
description of all the datasets (Model E: Table ^) . The ther- 
mal part of the hybrid electron distribution Comptonizes 
the disc blackbody photons into a broader soft component, 
while the non-thermal electrons Comptonize the disc pho- 
tons to make the harder spectrum which is reflected from the 
disc. The parameters of the thermal Comptonizing plasma 
are similar in all cases: kT e « 5 keV, r e = 1 — 2 (for a 
spherical geometry). They are consistent with the parame- 
ters derived above for a purely thermal origin for the soft 
excess (compPS or thComp), but are rather different from 
the parameters derived from fitting the high state spectrum 
of Cyg X-l, where kT c ~ 40 keV was required (Gierliriski et 
al. 1999). 

Another possibility to describe the complex soft compo- 
nent in these spectra is to assume a sum of a disc blackbody 
and an additional single blackbody. The resulting fits are 
good (model F: Table |2|), with the disc blackbody temper- 
ature at fcTdbb ~ 0.3 — 0.4 keV, with a higher temperature 
for the additional blackbody at fcTbb ~ 0.5 — 0.8 keV. 



Table 3. Parameters of the reprocessed component from fitting 
the double Comptonization model to GS 2000+25 data 

parameter Nov 5 Oct 18 Sep 8 

0.71±°;g o.52l° : °l 0.43i^ 8 

€ (18l 3 °) x 10 3 (10+ 25 4 ) x 10 3 (13+^) x 10 3 

Rm 15lf 5 8.2+ 3 'i 6+ 3 - 2 



The broadness of the soft component cannot be simply 
attributed to relativistic effects: the GRAD model, which in- 
cludes these effects, does not provide a good description of 
our data. The x 2 values shown in Table Q are usually even 
worse for the GRAD model (model G) than for the simple 
disc blackbody (model C). In these fits the inner disc radius 
has been fixed at the marginally stable orbit at 6R g . The fits 
can be improved if R[ n is a free parameters, in which case 
-Rin usually increases to > 20 7? g . The relativistic effects are 
not important then, and the model is equivalent to diskbb. 

The data clearly require a presence of the reprocessed 
component. To illustrate its significance we removed this 
component from the double Comptonization model (model 
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D), and fit parameters of the two continua, including kT c 
for the first component and Ah- This resulted in bad fits, 
as given as model in Table ^. The derived parameters 
for the reprocessed component are similar to those previ- 
ously inferred for soft states of black hole systems (Zycki 
et al. 1998; Gierlihski et al. 1999): the main features are 
that the reprocessor is highly ionized and it requires further 
smearing/broadening compared to predictions of the sim- 
ple model. This is demonstrated in Table ^| where we show 
parameters of the reprocessed component, as obtained from 
fitting model D to the Nov 5, Oct 18 and Sep 8 datasets. We 
model the smearing as due to relativistic and kinematic ef- 
fects, but spectral features can be additionally significantly 
broadened by Compton down- and upscattering, as the pho- 
tons diffuse through the ionized disc layers (e.g. Ross, Fabian 
& Young 1999). Therefore, R[ n as determined from our fits 
is likely to represent the lower limit on the inner radius of 
the reflecting disc. 

The bolometric luminosity of GS 2000+25 is rather 
low, if the distance to the source d = 2 ± 1 kpc (Calla- 
han et al. 1996) is adopted. Even for the earliest dataset 
analysed, Sep 8th, the bolometric luminosity, L — (1.1 — 
2.0) x 10 3r di erg s corresponds to less than 2 per cent 
of the Eddington luminosity. The luminosity on Dec 8th is 
~ 0.003 I/Edd, while one week later, during the first observa- 
tion in the LS the luminosity is ~ 2 x 10 -4 LecW- It is rather 
unclear whether this source was showing the same spectral 
characteristics as other sources (both transient and persis- 
tent), despite the low L/Z/Edd, or the distance is severely 
underestimated. 



3.2 GS 1124-68 

Ginga spectral data of GS 1124-68 (Nova Muscae 1991) have 
been analysed by a number of authors (Ebisawa 1991; Ebi- 
sawa et al. 1994; Esin et al. 1997). We have already anal- 
ysed a number of datasets in Zycki et al. (1998). In par- 
ticular, on 18-May-1991 (~ 150 days after the peak of the 
outburst) the source was observed in a soft state for the 
last time. Fitting the wabs(diskbb + thComp + rel- 
repr) model to the 1-20 keV data we obtain a very bad 
fit, Xv ~ 241/24, with overall absorption column Nu = 0. 
Replacing the diskbb model by Comptonized blackbody, 
modelled by a second thComp component, yields dramatic 
improvement, xi = 14.8/24, with TVh fixed at the inter- 
stellar value 1.6 x 10 21 cm -2 (note that in Zycki et al. 
1998 we used a power law to approximate the Comptonized 
hard component). Again the Comptonized soft component 
is rather steep, so although its slope is well constrained 
(r so f t = 5.6 ± 0.2), kT e and r e cannot be separately deter- 
mined in the approximate thComp models. Using CompPS 
for the soft component again demonstrates that the data re- 
quire that the Comptonizing medium has low temperature 
and high optical depth. 

The results are robust with respect to the model for the 
hard component: the wabs(diskbb + compPS + rel-refl 
model gives xi ~ 320/23, i.e. clearly an unacceptable fit. 
Replacing the diskbb component with thComp we again 
obtain a very good fit, 15.5/23. Thus the presence of the 
additional Comptonization on the soft component does not 
seem to be an artifact of the approximate character of the 
Comptonization models. 
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Figure 2. Upper panel: Unfolded spectrum of XTE J1550-564 
on Nov 22, 1998, during the long minimum after the peak of its 
outburst, showing Comptonized soft component, hard power law 
tail and X-ray reprocessed component. Lower panel: contributions 
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The hybrid Comptonization model gives a good fit to 
the data, with \ v = 17.1/25. The Maxwellian temperature 
is kT c = 5.40 ± 0.05 keV, the cloud optical depth is 2.94 ± 
0.10, similar to those derived for GS 2000+25. The model 
comprising a disc blackbody with an additional blackbody 
gives xt ~ 18.5/23. The model with the GRAD component 
for the soft emission does not give a good fit, xl = 215/25. 



3.3 XTE J1550-564 

The outburst of XTE J1550-564 in September 1998 and sub- 
sequent evolution were very well covered by pointed RXTE 
observations. Basic spectral analyses of the large collected 
dataset were performed by Sobczak et al. (2000 and refer- 
ences therein), while correlations between timing and spec- 
tral properties were studied by Homan et al. (2001). Wilson 
& Done (2001) studied in the source spectral evolution on 
the rising phase of the outburst, as the source made a transi- 
tion from the low state to the very high state. However, here 
we are interested in Intermediate State. We identify this as 
occurring during the long minimum after the peak and be- 
fore the source flux recovered again (see ASM lightcurve in 
fig. 1 in Homan et al. 2001), and we analyse a spectrum taken 
in the middle of this interval, on 22 Nov 1998 (MJD=51139). 

We fit the 2.5-30 keV data with the WABS* (diskbb + 
thComp + rel-repr) model, i.e. the hard component is 
modelled as a Comptonized emission with a corresponding 
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reprocessed component. When Nh is fixed at the likely inter- 
stellar value of 4 x fO 21 cm" 2 (Wilson & Done 2001), the fit 
is poor, xl = 79/54, with the strongest residuals between 5 
and 9 keV. Allowing for free Nn improves the fit somewhat; 
the best fit has xl = 72/53 for TVh = 0. Substituting now 
a second thComp component in place of the diskbb, yields 
a significant improvement, the best fit has xl ~ 41.7/53 
with Nn fixed at 4 x 10 21 cm" 2 (Fig. J). The hard compo- 
nent has a slope of ri mr d — 2.12l[5' < ( ) ) ;J and its reprocessed 
spectrum is detected with amplitude Q/2ir = I-IZiTm fr° m 
highly ionized £ = (30^23) X 10 4 material which is strongly 
smeared, corresponding to Ri n = 6 +3 R g (assuming the disc 
inclination of 60°). Removing the reprocessed component 
altogether gives a very poor fit with xl = 185/56. 

The parameters of the thComp model for the soft com- 
ponent are similar to those derived for GS 2000+25 and 
GS 1124-68 in that they give a steep power law tail to 
the disc blackbody (r so ft = 5.4 ± 0.4). This broadening 
of the soft component can be equally well matched by hy- 
brid Comptonization models (xl = 46.7/54) or as the sum 
of disc blackbody with an additional (hotter) blackbody 
(xl = 28.2/52). The fit with the GRAD model is rather worse, 
xl = 69.5/54 

Adding HEXTE data to the fit does not change our 
conclusion about the Comptonization of the soft compo- 
nent. With two COMpPS components we obtain a good fit, 
xl = 57/87, with the parameters of the soft component very 
similar to the previous case of PCA data alone. 

The distance to the source is rather uncertain, as 
is the mass estimate of the compact object and indeed 
its very nature. Estimates of the L/L-^dd ratio are there- 
fore very uncertain. The 2-30 keV flux in these data is 
~ 5.5 x 10 -9 erg cm -2 s _1 which is fa 30 times lower than 
that on 19th Sept (MJD = 51075), during the bright flare. 
However, our spectral analysis indicates that significant flux 
was emitted below 2 keV. The bolometric correction is rather 
uncertain, depending primarily on whether the seed photons 
are assumed to be simple blackbody or disc blackbody. We 
derive the bolometric (unabsorbed) , 0.1-1000 keV flux to be 
(1.4-2.7) x 10~ 8 erg cm -2 s" 1 w (3-5) x F(2-30keV). On 
the other hand, the bolometric correction for the flare spec- 
trum seems to be < 2, so the ratio Laaie/L ~ 10, possibly 
less. 

Adopting the distance d fa 2.5 kpc, based on optical 
observations during the 1998 outburst (Sanchez-Fernandez 
et al. 1999), we obtain rather low peak luminosity, Lfl arc fa 
2xl0 38 erg cm" 2 s~\ hence L fa 2xl0 37 erg cm" 2 s" 1 . The 
latter corresponds to L/Lsdd ~ 0.03 df.s-Mg" (where Me = 
M/(6Mq)), similar to the soft-hard transition luminosity 
for Cyg X-l (Gierliriski et al. 1999; Di Salvo et al. 2001). 



4 THE HIGH AND VERY HIGH STATE DATA 

A number of SXTs showed spectral/timing states classified 
as a VHS. Nova Muscae 1991 showed such a behaviour dur- 
ing the first ~month around the peak of the outburst (16th 
Jan 19991) (Ebisawa et al. 1994; Takizawa et al. 1997; Bel- 
loni et al. 1997; Rutledge et al. 1999). GRO J1655-40 went 
through VHS during its outburst in March - September 1997 
(Mendez, Belloni & van der Kils 1997). XTE J1550-564 dis- 
played a variety of spectral/timing states including the VHS 



(Homan et al. 2001; Wilson & Done 2001). A peculiar spec- 
trum was shown by GS 2023+338 at the peak of its 1989 out- 
burst: the spectrum was very hard, with high energy cutoff 
at the low value of ~ 20 keV (Zycki, Done & Smith 1999b). 

4.1 GS 1124-68 

4.1.1 VHS - Jan 11th data 

During this observation the source was still on the ris- 
ing phase of the outburst, but its bolometric luminosity 
(with the best model described below) was fairly high, 
L fa 0.25I/Edd rf|.5 M^ 1 (F = 1.4x 10~ 7 erg cm" 2 s" 1 ). Sim- 
ilarly to the case of the IS spectra, we find here that the 
soft component cannot be described by the disc blackbody 
model. More specifically, the best fit of the wabs (diskbb + 
thComp + rel-repr) model has xl = 184/31. Replacing 
the disc blackbody component by the Comptonized black- 
body model thComp, gives us a good fit, xl = 27.7/30. 
This confirms our earlier result (Zycki et al. 1999b, see cap- 
tion to fig. 5), that the soft component in that spectrum 
is Comptonized, in addition to the obvious Comptonization 
of the hard component (note that our previous result was 
obtained with a power law model for the hard component, 
while here we use a Comptonized spectrum for it). The spec- 
trum is plotted in Fig. |§| Reflection of the hard compo- 
nent is significantly present in the spectrum, its amplitude 
is Q./2tt = 0.24+°' 2 3, it is strongly ionized £ = (3t$. 7 ) x 10 3 
and additionally smeared. The smearing corresponds to in- 
ner disc radius of Ri n = 16l|. The model with no repro- 
cessed component gives xl = 555/33. The spectral slope of 
the soft Comptonized component is F so f t = 2.94lQ 2g, rather 
flatter than seen in the IS spectra above. With thComp or 
COMPTT, we can only give the lower limit to the plasma 
temperature, kT c > 2 keV. The soft component can also be 
well described by the sum of disc blackbody (fcTdbb ~ 0.5 
keV) and and additional, hotter (fcTbb ~ 1.2 keV) black- 
body, giving the overall xl ~ 13.0/31. We were not able 
to find a good fit with the hybrid Comptonization model 
(xl > 150/32). Similarly, the GRAD model used for the soft 
component gives the overall poor fit with xl = 239/33. 

Note that the overall spectrum is rather hard (F ~ 2) 
and it lacks a dominant soft component, making it rather 
similar to typical low state spectra. This led Miyamoto et 
al. (1995) to suggest that the transition from spectral hard 
to soft state occurs at much higher luminosity on the rising 
phase of the outburst than on the declining phase, i.e. the 
source shows a kind of hysteresis. However, with more de- 
tailed spectral fitting it is evident that the VHS spectrum is 
rather different from typical low/hard state spectra: the soft 
component is not a pure disc blackbody emission and the re- 
processed features correspond to highly ionized plasma. On 
the other hand, during the decline, spectra at similar lumi- 
nosity level were totaly dominated by the soft component 
(see discussion below and Fig. |(j and the Jan 11. vs. Feb 
13. spectrum). Thus indeed the developement/destruction 
of the optically thick disc spectra does not seem to be sym- 
metric in m, in a time-dependent situation. 
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Figure 3. Spectra of GS 1124-68 (Nova Muscae 1991) in Very High and High states. The thick solid histogram represents the best fit 
total model, the histograms drawn in thin line plot the model components: solid histogram - the soft component, dashed histogram 
— the hard component, dotted histogram — the reprocessed component. In all shown spectra the soft component is significantly better 
described by a Comptonized (disc) blackbody than the simple disc blackbody emission. On Jan 11 the source was still on the rising 
phase of the outburst. The source luminosity peaked on Jan 16 (the hybrid Comptonization model is plotted for these data). On Feb 13 
the bolometric luminosity is very similar to that on Jan 11, but the two spectra is rather different (Miyamoto et al. 1995). In the last 
spectrum shown the hard component is negligible. 



4.1.2 VHS - Jan 16th data 

On Jan 16th the source's flux reached maximum. The spec- 
trum is dominated by a soft thermal component with a weak, 
steep power law tail (Ebisawa et al. 1994). 

Fitting the wabs (diskbb + thComp + rel-repr) 
model to 1-30 keV data we obtain a poor fit, xt = 73.8/32, 
for JVh fixed at the interstellar value 1.6 x 10 21 cm -2 . Allow- 
ing JVh to be free does not improve the fit. The strongest 
residuals are below ~ 5 keV, again indicating a more compli- 
cated shape of the soft component. Replacing the disc black- 
body with a Comptonized component gives an extremely 
good fit: wabs (thComp + thComp + rel-repr) (with 
blackbody seed photons) has xt = 11-7/32 (with Nn fixed). 
However, as opposed to the Intermediate State spectra, here 
the curvature of the soft component is clearly seen and, 
consequently, the electron temperature is very tightly con- 
strained at kT c = 1.12±g;{£ keV. The Thomson depth of the 
Comptonizing plasma is large, with r c ~ 8 for slab geometry. 

Again, an acceptable fit can also be obtained by a hy- 
brid electron distribution: diskbb+ compPS + rel-repr 
gives xt ~ 28.3/32 assuming that the seed photons for 
the compPS are the disc blackbody, although this requires 
JVh larger than the interstellar value, Nn = (3.5 ± 0.8) x 
10 21 cm -2 . The thermal part of the hybrid electron distri- 
bution is then again a rather cool, kT e = 6 keV, optically 
thick plasma, r c = 2.3. Alternatively the soft component 
can be modelled as a sum of a simple blackbody and a 
disc blackbody. This gives xi = 12.1/30 and parameters: 
fcTbb ~ 0.50 keV and kT dbb fa 1.1 keV, thus the blackbody 
seems to have a lower temperature than the disc blackbody 
component, opposite to what was observed for IS spectra. 



Energetically, the blackbody component contributes ~ 30 
per cent of the total flux of the soft component. The GRAD 
model for the soft component does not provide a good fit, 
giving the overall xl = 82.4/33. 

The dominance of the soft component in these spec- 
tra give perhaps the best view of its shape. Hence we have 
also fit the data a model of multi-temperature blackbody 
emission with the radial dependence of the temperature dif- 
ferent from the standard formula, T(r) oc r -3 ^ 4 . Our moti- 
vation for such a model comes from recent work by Watarai 
et al. (2000), who draw attention to the fact that in lu- 
minous accretion discs, M ~ MEdd, the T(r) dependence 
is flatter than the standard one due to a fraction of en- 
ergy being advected radially. We have tried both a simple 
model, with T(r) = Ti n (r/Rm)~ a (Tin and a free param- 
eters), and a more complex one, where T(r) was actually 
computed from the radial structure of a vertically averaged 
Shakura-Sunyaev (1973) disc, with the radial advection in- 
cluded. However, neither of the two models (together with 
thComp + rel-repr for the hard component) provides a 
good description of the data. The broadening detected in the 
soft component is much more evident that can be produced 
by such models. While it is to be expected that the simple 
disc blackbody approximation breaks down when the source 
luminosity reaches ~ L^dd (see, e.g., Beloborodov 1998), 
simple 'first order' corrections to the disc blackbody models 
are not sufficient. 

Note that in Zycki et al. (1999b) we have analysed 
this dataset and concluded that it is compatible with pure 
disc blackbody emission with a hard tail. i.e. no additional 
Comptonization of the soft component was required to make 
the fit statistically acceptable. That conclusion resulted, 
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however, from our using larger systematic errors of 1 per 
cent. More specifically, assuming the 1 per cent system- 
atic error, the model diskbb+ thComp + rel-repr gives 
xl = 33.0/32. However, the model thComp + thComp + 
rel-repr gives xl = 7.5/30 when fit to the same data. The 
conclusion that the Comptonized blackbody model gives a 
better fit than the disc blackbody is obviously independent 
of the assumed value of systematic errors. 

4.1.3 VHS - Feb 13th data 

This is also a Very High state spectrum, at a luminosity 
lower by a factor of ~ 2 than the peak luminosity. Here we 
also find evidence for a Comptonized soft component. 

The model wabs(diskbb + thComp + rel-repr) 
gives xl = 119.3/25 for fixed Ah- Letting Ah free gives 
xl = 63.8/24 and N H = 0+ 03 x 10 21 cm~ 2 , clearly less 
than the interstellar value. Replacing diskbb with the rel- 
ativistic model GRAD improves the fit somewhat, giving 
xl = 38.8/24. 

The Comptonized blackbody model with a hard tail 
gives a very good fit, xl = 8.0/23, using the thComp model 
for the soft component, or xl = 9.8/23 using the COMPTT 
model. Again the data require the additional Comptonizing 
medium to be at fairly low temperature kT c = 0.851q'q5 
keV, and optically thick r c = 9.4 ± 1.3 (for slab geometry). 
Reprocessed emission from the hard component is required 
in the spectrum, with amplitude Q,/2tt — 0.51q'| and ioniza- 
tion parameter £ = (2.3^21) x 10 4 - The relativistic smear- 
ing is not required by the data, but they are consistent with 
smearing with _Rj n > 17i? g . Setting the reflection amplitude 
S1/27T = and fitting all the other parameters of the contin- 
uum gives a bad fit with xl = 95/26. 

Again, an alternative description of the soft spectral 
broadening can be found using diskbb with an additional 
blackbody {xl = 16.1/22). However, the hybrid Comp- 
tonization does not provide a good fit. 

4-1-4 High State - March 18th data 

Approximately two months after the peak of the outburst 
Nova Muscae 1991 made a transition to the High State, 
where its spectrum showed a strong soft component with 
only very weak hard tail, if at all present (Ebisawa et al. 
1994; Esin et al. 1997; see Fig. [|. Ebisawa et al. (1994) no- 
ticed that the disc blackbody description failed completely 
for these data (xl ~ 10). We have modelled the spectrum 
obtained on 18th March. Applying the diskbb + thComp 
model gives xl = 54/26. The fit is bad, since the Comp- 
tonized component is steep, F ~ 8, so it makes a significant 
contribution to the soft component and, as a consequence, 
misses the points at E > 10 keV. Thus again the soft com- 
ponent has to be described by a model broader than the 
multi-colour disc blackbody. 

This can also be demonstrated by ignoring data above 
10 keV, thus leaving only the soft component (contribution 
from the hard component to E < 10 keV is negligible) and 
attempting to model it by a disc blackbody spectrum. The 
fit is very bad, xl = 756/12, even if A'h is allowed to vary. 
The grad models gives a better fit, xl = 329/12, but this 
is still far from satisfactory. These data can be well fit by a 
Comptonized disc blackbody, giving xl = 10.9/10. 



The overall spectrum can be well fit using the COMpPS 
model with a hybrid electron energy distribution, xl = 
18.0/24. The low energy, thermal (kT c « 3 keV) electrons 
Comptonize the seed, disc blackbody photons broadening 
somewhat the soft component, while the power law tail in 
the electron energy distribution gives rise to the weak power 
law tail seen above 10 keV. An additional blackbody also 
gives a good fit: wabs (diskbb + bbody + thComp) gives 
Xl = 14.6/26. 

4.2 Very High State of GRO J1655-40 

GRO J 1655-40 went through a sequence of spec- 
tral/temporal states during decline after the outburst ob- 
served in March - September 1997. Analysis of the data by 
Mendez et al. (1997) suggests that during the May 28 ob- 
servation the source was in the VHS: the energy spectrum 
contained both a strong soft component and a power law 
tail. Likewise, the PSD appears to be a sum of a power law 
component and a flat-topped component (see Takizawa et 
al. 1997). 

We use the 2.5-30 keV top layer RXTE/PCA and 30- 
200 keV RXTE/HEXTE data for analysis. The N H value 
towards this object is rather uncertain: the E(B — V) = 1.15 
(Bailyn et al. 1995) corresponds to Ah,e(b-v) = 6.5 x 
10 21 cm -2 for standard extinction curve. Fitting the ASCA 
data of the August 1995 observation Zhang, Cui & Chen 
(1997a) found A'h = 8.9 x 10 21 cm~ 2 , using a power law 
model for the hard component, while Gierlihski, Maciolek- 
Niedzwiecki & Ebisawa (2001) found A'h ~ 7 x 10 21 cm" 2 
using a Comptonization model for the hard component. Fit- 
ting the same data, we find A'h values both higher and lower 
than 6.5 x 10 21 cm -2 depending on the model, therefore we 
let this parameter be free in our fits. We assume the source 
inclination is i = 70°. 

We begin with PCA data alone. Unusually, the simple 
model wabs (diskbb + thComp + rel-repr) gives a good 
fit, xl = 36.2/54, with A'h = (4.3 ± 1.4) x 10 21 cm -2 . Fixing 
Ah at Ah,e(b-v) gives xl = 42.5/55. Further broadening 
of the soft component is allowed by the data although not 
required. Replacing the disc blackbody model with thComp 
yields X l = 33.0/53, again for a cool kT c ~ 1 keV, optically 
thick, r c ~ 10 (for a slab geometry) plasma, Comptonizing 
seed photons at a rather low temperature To ~ 0.3. These 
results are very similar to those of Zhang et al. (2000), who 
fitted similar models to the ASCA data of GRO J1655-40 
and GRS 1915+105. Comptonization of the soft component 
in this source, in high luminosity states, is also reported by 
Kubota, Makishima & Ebisawa (2001). 

Again, the alternative descriptions of the broad soft ex- 
cess give excellent fits to the data. A hybrid electron distri- 
bution gives xl = 33.2/55, while an additional blackbody 
gives xl — 35.2/53, both fits with A'h = Ah.e(b— v) ■ 

The central black hole in this object is postulated to 
be a rotating one based on the observed high temperature 
and small amplitude of the thermal disc emission (Zhang 
et al. 1997a). Gierlihski et al. (2001) infer black hole spin 
a — 0.7 — 0.8, fitting a relativistic disc spectrum model, with 
photon propagation in the Kerr metric to the 1995 ASCA 
data. However, in order to obtain a good fit, Gierlihski et al. 
need to introduce a smeared edge at 2-3 keV, to model dis- 
tortions present in the spectrum. Replacing the disc black- 
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Figure 4. Model fit to the VHS RXTE data of GRO J1655-40 
( x 2 = 90/113). The model comprises the hybrid Comptonization 
and additional thermal Comptonization. The reprocessed compo- 
nent is clearly detected, with amplitude Q/(2ir) fa 1, high ioniza- 
tion, £ sa 10 s and relativistic smearing, _Rj n as 16R g . 



body model in our fits with a disc blackbody in Kerr metric 
we obtain xt = 44.6/54, fixed at -/V h ,e(b-v), but with the 
inner disc radius and central mass as free parameters (to 
compensate for the fixed value of the black hole spin in our 
model; see discussion in Gierlihski et al. 2001). Thus using 
the Kerr metric disc model does not improve the fit com- 
pared to the simple disc blackbody model. In view of the 
complex soft component shape which is generally seen in the 
soft spectral states we caution strongly against interpreting 
any broadening of the thermal emission solely as relativistic 
effects (although such effects should be present). 

The above models are no longer adequate when the 
HEXTE data are added. An additional, harder continuum 
component is required to model the broadband spectrum, 
as already found for the VHS spectra of XTE J 1550-564 by 
Wilson & Done (2001). A possible description of the spec- 
trum is provided by the hybrid Comptonization model with 
an additional, thermal Comptonization. The slopes of the 
two Comptonized tails (the non-thermal tail and the addi- 
tional Comptonization) are not well constrained individu- 
ally, so we fit the model fixing the slope of the non-thermal 
electron distribution at 7 = 5 (i.e. steeper than in the fit to 
PC A data alone). We allowed both continua to give rise to 
the reprocessed component, but with the same parameters 
(Q/(2n), £ and Ri n ). The overall model, WABS(cOMPPS hy brid 
+ rel-refl + thComp + rel-repr), gives a good fit to 



the data, xt = 90/113 (Fig. ^|). The reprocessed component 
is clearly detected in the spectrum: without it the fit is much 
worse, xi ~ 158/116. The reflection amplitude Q/(2tt) ~ 1 
(lower limit ~ 0.5 but with poorly constrained upper limit), 
the ionization parameter xi ~ 10 , and the inner disc radius 
-Rinl6i? g . This differs from the results of Gierlihski et al. 
(2001), who did not find any significant reflection features 
in ASCA data of this source. 



5 DISCUSSION 

We have re-analysed archival data of a number of SXT 
sources in soft spectral states in an attempt to character- 
ize their spectral properties in more detail than has gener- 
ally been done so far. In most previous analyses of IS, HS 
and VHS spectra the assumed spectral model consisted of a 
disc blackbody emission, a power law, a smeared edge and 
a gaussian line near 5-9 keV. Of these, the edge and the 
line are the least accurate representation of the underlying 
physical processes. We replace them with a more physical 
model of X-ray reprocessing by an optically thick plasma. 
Also, a power law is an inaccurate approximation to Comp- 
tonized emission at energies close to the seed photons as true 
Comptonized spectra have a low energy cutoff. In the soft 
spectral states the seed photons from the disc are within the 
observed bandpass, so such effects become important. 

Two results from the above data analyses seem to be 
robust: firstly, the soft component cannot be adequately 
described by the simple disc blackbody model and, sec- 
ondly, the reprocessed component is present in the data, 
it is strongly ionized and smeared. 

The complexity of the soft component is not merely an 
artifact of using the approximate diskbb model. Including 
the torque- free boundary condition, colour temperature and 
relativistic effects in photon propagation does not substan- 
tially change the overall shape of the spectrum, although 
the derived model parameters can be significantly different 
(e.g. Ebisawa 1991; Gierlihski et al. 1999; Merloni, Fabian 

6 Ross 1999). Direct fitting of our relativistic disc model 
GRAD, which includes the corrections, supports our claim. 
Although the best fits are somewhat better than the model 
employing diskbb, they are much worse than fits employ- 
ing Comptonized black body. Similarly, using the Kerr met- 
ric disc emission model did not provide better description 
of the spectrum of GRO J1655-40 compared to the simple 
disc blackbody model. The soft component is substantially 
broader than optically thick disc models, irrespective of rel- 
ativistic corrections. 

Using the Gmga or RXTE data we are not able to 
differentiate between various possible modifications of the 
disc blackbody spectrum which ensure a good description 
of data. The two main possibilities are additional Comp- 
tonization of the disc emission or an additional blackbody 
component. We discuss each of these in turn. 



5.1 The complex soft component: additional 
blackbody 

One possibility to explain the complex shape of the soft com- 
ponent is to assume that it is a sum of a disc blackbody and 
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additional blackbody emission, indicating that some frac- 
tion of the disc material is heated to rather higher temper- 
atures. Physically we might expect this if the hard X-ray 
power law arises from magnetic reconnection. The regions 
of the disc below the intense flares are then heated by the 
irradiating flux and the thermalized (non-reflected) radia- 
tion would be emitted as a quasi-blackbody spectrum at 
temperatures somewhat higher than the mean disc emission 
(Haardt, Maraschi & Ghisellini 1994). Alternatively, for a 
central source geometry, the inner disc is additionally heated 
by irradiation. In all these models the additional blackbody 
is powered ultimately by the hard component. 

Energetically this is certainly feasible for the IS spectra. 
Consider, for example, the GS 1124-68 May 18th data. The 
bolometric (0.01-1000 keV) flux from this model is (in units 
of 10 -9 erg cm -2 s _1 ): F to t w 25, the hard X-ray flux is 
^hard ~ 2.0, while the flux in the additional blackbody is 
Fbb = 0.88 w 0.44Fi mr d. The reprocessor albedo is observed 
to be a ss 0.75 (ionization of £ m 3000, so that FeXXV and 
FeXXVI are the dominant ionization stages of iron). Thus 
the thermalized fraction is expected to be 1 — a m 0.25. Given 
the uncertainties in the modeling this is probably consistent 
with the observed fraction « 0.4. 

The observed flux F hh ^ 9 x 10~ 10 erg cm -2 s _1 
at Tbb ~ 0.7 keV implies an emitting area S = 
-iird 2 F^/(aT 4 cosi), where a is the Stefan-Boltzmann con- 
stant. This gives a tiny area, S w 10 13 cm 2 rss 13 Rg, assum- 
ing distance d = 3.5 kpc, M = 6M0 and neglecting the 
colour temperature correction. This is consistent with mag- 
netic flare illumination, if the illuminating area is small, i.e. 
if the flares are located not too high above the disc (Nayak- 
shin & Kazanas 2001). In the central source geometry the 
distribution of irradiation flux is somewhat steeper than r~ 3 
close to the inner disc radius (see fig. A2 in Zycki et al. 
1999a). This gives an excess of emission, if modelled by the 
diskbb model, but the excess is not sufficient to account 
for the observed shape of the soft component. Partial over- 
lap between the central source and the disc might solve the 
problem. 

However, the soft component is overwhelmingly domi- 
nant in the HS spectra and reprocessing of the very small 
hard X-ray flux is not consistent with the flux required for 
the additional blackbody. An extreme example here is the 
HS GS 1124-68 March 18th dataset, where the hard X-ray 
(E > 10 keV) flux is less than 1 per cent of the total flux, 
while the additional blackbody component contains ~ 1/3 
of the total flux. The situation is not so extreme for the 
VHS spectra, where there is some hard X-ray emission. The 
observed fluxes of hard X-rays are rather smaller than the 
fluxes of the additional blackbody components, which rules 
out the possibility of the blackbody being the thermalized 
hard X-ray illumination, if isotropic emission of the latter 
is a ssume d. The peak spectrum of GS 1124-68 (Jan 16th; 
Sec 4.1.2) is a good example here: the fluxes in the 0.01-10 3 
keV band (in units of 10 -8 ergs cm -2 s _1 ) are: F to t = 28.3, 
F disc bb = 17.1, Fhard = 3.4 and F bb = 7.8 > F hard . If the 
hard X-rays were beamed towards the disc, both the ther- 
malized emission and the reprocessed component would be 
enhanced. The amplitude Q./2n is not always well constained 
in VHS data due to complex continuum curvature; for ex- 
ample for the Jan 16th data Q./2n = 1 
Feb 13th = 0.5±°*. 



5.2 The complex soft component: Comptonization 

The two-component description of the disc emission does not 
seem appropriate for the brighter spectra, especially those 
lacking a strong hard component. The additional Comp- 
tonization of the disc emission provides an alternative de- 
scription to those (and all other) spectra. The parameters 
of the Comptonizing plasma are not well constrained, but 
generally a cool, kT c = 2 — 4 keV, optically thick, r c ~ 5, 
plasma gives the best fit. The main difference between this 
and the same spectral index Comptonization from a hotter, 
optically thinner plasma is in the possible contribution of 
the seed photons to the final spectrum for r e < 1. However, 
this subtle difference is not sufficient to uniquely select the 
best fit solution, given the rather poor energy resolution and 
low energy bandpass of the data used here. 

We can envision a number possibilities for the origin of 
the Comptonization of disc emission. 



5.2.1 Warm atmosphere 0} the accretion disc 

The soft photons from the disc interior can be Comptonized 
as they escape through a hot, ionized topmost layer of the 
disc. 

If the Comptonizing layer is heated by external irradi- 
ation only, i.e. viscous heating is negligible, its maximum 
temperature is the inverse Compton temperature, Tic, at 
which Compton heating is balanced by Compton cooling. 
This temperature can be estimated from 



4fcTi C / ue dE = / Eu E o- K N(E)dE, 



±f%, while for the 



(1) 



where ue is the radiation energy density, and ctkn(-B) is the 
Klein-Nishina cross section. For a spectrum consisting of a 
strong soft component and a steep power law, Tic is not 
much larger than the temperature of the soft component. 
For example, for the May 18th data of GS 1124-68 we calcu- 
late Tic ~ 0.9 keV, and similar values are obtained for the 
other IS spectra. This is lower than the inferred temperature 
of the Comptonizing plasmas, implying that there must be 
additional heating of the plasma through e.g. viscous dissi- 
pation (Janiuk, Czerny & Madejski, 2001). Only for the 11th 
Jan spectrum of GS 1124-68 is Tic ~ 4 keV comparable to 
the required kT c of the Comptonizing plasma. 

This picture then implies a three-phased accretion sce- 
nario: the hot, kT c ~ 100 keV, plasma with the Comptoniza- 
tion parameter y ~ 1, producing the hard X-ray continuum; 
the warm, kT e — 1-10 keV plasma with y ~ 0.3, responsi- 
ble for the additional Comptonization of the disc emission; 
and the optically thick, 'cold', kT c ~ 0.1-1 keV disc. A sim- 
ilar accretion scenario was also suggested for the Seyfert 1 
galaxy NGC 5548 (Magdziarz et al. 1998). The clear soft 
X-ray excess component in that object could be well mod- 
elled by a Comptonization of accretion disc photons in a 
warm (kT c ~ 300 eV), optically thick (r c ~ 30) plasma. 
Magdziarz et al. speculated that the plasma could form a 
transition region between the outer, cold accretion disc and 
the inner hot flow. Similar Comptonized soft X-ray excesses 
are present in two Narrow Line Seyfert 1 galaxies: PG 1211- 
143 (Janiuk et al. 2001), and PKS 0558-504 (O'Brien et al. 
2001). A possible analogy of such three-phased configura- 
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tion around accreting black holes with the Solar corona was 
suggested by Zhang et al. (2000). 

However, the warm plasma cannot form a continuous 
later covering the cold disc. It is optically thick and com- 
pletely ionized so would produce no spectral features, con- 
trary to results of our data analyses. More robust (given 
the uncertainties in optical depth and temperature) argu- 
ments against such a geometry are provided by considering 
the energy balance between the cool disc and warm skin. 
Following the approach of Haardt & Maraschi (1991), the 
energy balance can be written as 

l~f + Vf = Fseed, (2) 

while for the warm skin 

AF Bcc d — ^compt — (3) 

where / the fraction of energy dissipated in the warm skin 
(disc dissipation 1 — /), r\ is the fraction of Comptonized 
emission returning to the cold disc, A is the Comptonization 
amplification factor, A = F CO mpt/ -fsccd, and all the fluxes 
are normalized to the total gravitational energy dissipated 
in the two layers. Solving Eq. (Q) and Eq. (|j|) , we obtain the 
condition for the existence of solution, 

r\A < 1. (4) 

Curiously, the Comptonization parameters from the fits to 
GRO J1655-40 imply A « 1.24 and n ~ 0.82, thus \r)A - 
1| < 0.02. This, in turn, implies / very close to 1, i.e. the 
entire energy dissipation in the warm skin, with the cold 
disc merely reprocessing and thermalizing the intercepted 
radiation. 

The situation is analogous to that with a continuous 
hot corona overlapping an optically thick disc (Haardt & 
Maraschi 1991). Even in the extreme case of / = 1 the 
produced spectra are not harder than T ~ 2, because the 
geometry alone implies A < 2. Any harder spectra imply 
e.g. a patchy corona, coronal outflow, etc. (e.g. Stern et al. 
1995; Beloborodov 1999). 

A possible solution to the reprocessed features problem 
is to assume that the scattering layer appears only in the 
inner part of the disc, close to the radius of maximum ef- 
fective temperature. The reprocessing features could then 
be produced further out, where the temperature is lower en- 
abling iron to recombine to H- and He- like ions, as observed. 
The relativistic smearing of the iron features, usually corre- 
sponding to Ri n > 10 R g supports this possibility. 

5.2.2 Hybrid electron energy distribution in the 
Comptonizing plasma 

Another possibility for the Comptonization of the soft pho- 
tons is to assume a hybrid, thermal/non-thermal distribu- 
tion of electron energies in an accretion disc corona (see 
Coppi 1999 for extensive discussion of relevant physical pro- 
cesses). The harder power law spectrum extending above 
~ 10 keV would then be a result of single inverse Compton 
upscattering of the disc photons by energetic electrons from 
the power law tail of the distribution. The Comptonized tail 
of the soft component is a result of scattering of the pho- 
tons by the low energy electrons, with roughly Maxwellian 
distribution. Such model was successfully applied to the soft 
state data of Cyg X-l (Gierlinski et al. 1999), where the soft 



component shows a clear excess above a disc blackbody com- 
ponent, while the hard power law extends without a break 
to at least 1 MeV (Grove et al. 1998). Fitting the CompPS 
Comptonization model (Poutanen & Svensson 1996) with 
such hybrid electron distribution we found that indeed it 
provides a good description of some of the spectra: the IS 
spectra of all three analysed objects are well described by it, 
as well as some VHS spectra: the peak spectrum of GS 1124- 
68 and the spectrum of GRO J1655-40. On the other hand, 
we were not able to find an acceptable description of the Jan 
11 and Feb 13 of GS 1124-68 with this model. 



5.2.3 Weak thermal corona 

Yet another possibility for the Comptonization of soft pho- 
tons is a weak, accreting corona, similar to that studied by 
Witt, Czerny & Zycki (1997), and Janiuk, Zycki & Czerny 
(2000). In these models the disc-corona transition is gov- 
erned by the condition of thermal equilibrium between the 
two phases of accreting plasma. This results in a rather weak 
corona, dissipating up to ~ 40 per cent of gravitational en- 
ergy. Details of the solution depend upon the role of advec- 
tive energy transfer (compare solutions in Janiuk & Czerny 
2000 and Janiuk et al. (2000), but generally the corona is 
optically thin, r c < 0.15 and its temperature varies between 
« 70 keV at R w 10 R e to 150 keV further away. The result- 
ing spectra are rather steep and roughly match the Comp- 
tonized excesses in the IS spectra, but such a weak corona 
would not be able to account for the strong Comptoniza- 
tion observed in e.g. Jan 16th spectrum of GS 1124-68. The 
model cannot account for the origin of the hard component 
seen in IS and VHS. There is, however, growing evidence 
of a very close link between hard X-ray and radio emission 
(see Fender 2000, 2001 for reviews). If the two are indeed 
linked, one may envision the hard X-rays originating from 
e.g. a base of a jet. 



5.3 X— ray reprocessing 

Spectral features near 6-9 keV are commonly observed in 
spectra of accreting compact objects, and they are indeed 
present in the spectra analysed in this paper. Reprocess- 
ing of hard X-rays by an optically thick plasma (Compton 
reflection, photo-absorption, fluorescent and recombination 
emission) is the standard interpretation of these features. It 
has been tested on numerous spectra of all kind of accret- 
ing objects, although most of the modeling work has been 
done with low/hard state spectra. The problem with soft 
states spectra is the complicated shape of the continuum 
emission near 5-10 keV, where the soft and hard component 
contribute roughly equally. In our work we do see that the 
parameters of the reprocessed component in the soft states 
spectra are dependent on the continuum model adopted. 
However, two results are robust from the work done so far: 

• the very presence of the reprocessed component, i.e. we 
have not been able to find a continuum description that 
would make the reprocessed component unnecessary, 

• high level of ionization of the reprocessor: again, irre- 
spectively of the continuum model, the reprocessing features 
correspond to highly ionized plasma, where most abundant 
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Figure 5. Effective absorption cross-section (free-free process 
only) as a function of energy, for a number of values of the gas 
temperature, T, and H number density, njj. Values of T and 
log(nH) are plotted next to corresponding curves. Since the ratio 
r cff/ T cs < 1 for E > 1 keV, the incident X-rays will not be ther- 
malized in the upper disc layers, where the process of Compton 
reflection is taking place. 

Fe ions are more strongly ionized than FeXVI (M-shell elec- 
trons removed), and in most cases the He- and H-like ions. 

The third result is that the spectral features are further 
broadened and smeared compared to predictions of the sim- 
ple model. We model the smearing as the relativistic Doppler 
effect and gravitational redshift, but it has to be kept in 
mind that Comptonization in the reprocessor is quite signif- 
icant, especially for strong ionization (e.g. Ross et al. 1999). 
Therefore the quantitative determinations of the inner ra- 
dius of the reprocessing accretion disc may be overestimated. 

A potential alternative to the reprocessing scenario 
might be absorption/emission in winds which are almost 
certainly launched from the vicinity of central objects ac- 
creting at rh close to the Eddington limit (e.g. Proga 2000 
and references therein). The characteristic 'P-Cygni' emis- 
sion/absorption profiles might turn out to be compatible 
with the data. This idea has not been explored in any de- 
tails and its quantitative tests are not possible at present, 
as the spectral predictions of the models are not ready yet 
to be quantitatively confronted with the data. 

One potential problem that might affect our conclusions 
is that the computations of the reflected continuum used in 
the rel-repr model (Magdziarz & Zdziarski 1995) assume 
that the only X-ray absorption process is the photo-electric 
absorption. However, the density in accretion discs around 
stellar mass black holes is high enough for the free-free pro- 
cess to become important. Another way of formulating the 
problem is to notice that - by the very definition of a soft 
state - the X-rays at E ~ 0.5 keV are thermalized. If so, 
then the reflected continuum could have a rather different 
shape at low energies than that used in our model. Precise 
computations of the reflected continuum in such conditions 
are not possible at present, as no photo-ionization code can 
handle densities ~ 10 20 cm -3 and thermalization of incident 
X-rays. We can however estimate whether the incident X- 
rays at E > 1 keV are likely to be thermalized or not. Solv- 
ing the vertical structure of an aP disc around 10 Mq black 
hole at R = 10 R g (see e.g. Rozanska et al. 1999) we find the 
H number density at the disc equator riH ~ 10 22 cm -3 . We 



compute the free-free absorption cross-section (Rybicki & 
Lightman 1979) as a function of energy for a number of val- 
ues of gas temperature and density. The results, plotted as 
the effective absorption, « c ff = W K a bs(K a bs + k C s) in units 
of the Thomson cross section as a function of energy, are 
shown in Fig. [s| Since the effective absorption coefficient is 
smaller than 1 for E > 1 keV, the incident X-rays will not 
be thermalized in the uppermost r cs ~ 1 layer, where Comp- 
ton reflection takes place. Thus the adopted model for the 
Compton reflected continuum is appropriate. 

5.4 Disc stability 

Perhaps the most puzzling behaviour of SXTs, in the con- 
text of the a-discs model (Shakura & Sunyaev 1973), is the 
strong soft thermal emission observed when the sources ra- 
diate at a significant fraction of the Eddington luminosity, 
when the aPtot-discs are unstable (Kato et al. 1998), and 
so perhaps they should not exist. On the other hand, these 
discs are stable for low L/Z/Edd (LS) when, observationally, 
they seem not to be present in the central regions of GBH. 
While the latter problem seems to have found its solution 
in the mechanism of disc evaporation (Rozanska & Czerny 
2000 and references therein), the former one is far from any 
satisfactory explanation. 

Here we construct the rh— S (accretion rate - surface 
density) diagram and mark the position of GS 1124-68, to 
quantitatively demonstrate the extent of the instability. This 
extends studies done by Gierlihski et al. (1999), who argued 
that in the soft state of Cyg X-l the accretion disc could 
remain on the stable branch, down to the last stable or- 
bit. However, we note that the soft state of Cyg X-l seems 
to correspond to an Intermediate State of GS 2000+25 or 
GS 1124-68 rather than a true High State, since its bolomet- 
ric luminosity is rather low, only « 0.03 LEdd (Gierlihski et 
al. 1999) , not much higher than a typical low state luminos- 
ity (Zhang et al. 1997b; Di Salvo et al. 2001). We estimated 
the bolometric luminosities of GS 1124-68 from models us- 
ing either disc blackbody or Comptonized disc blackbody 
spectra, which should give an upper limit to the bolometric 
correction. We further assumed a 6 Mq black hole at d — 3.5 
kpc. Since the largest value of d found in literature is ~ 5 
kpc (Esin et al. 1997), the uncertainty of d 2 should be less 
than factor of 2. Figure ^ shows positions of the source on 
the m-E diagram, demonstrating that the disc was stable 
in LS, but it was unstable (at least in some inner region) 
whenever the source was observed in a soft state. One might 
then speculate that the result of the instability is a slight 
reconfiguration of the disc, so that the disc emission is still 
optically thick, but with the additional Comptonization in 
some warm plasma. There is a number of instabilities operat- 
ing in radiation pressure dominated discs (Kato et al. 1998; 
Gammie 1998; Blaes & Socrates 2000), and it is possible 
that their outcome may be a clumpy accretion flow (Krolik 
1998). On the other hand, the amplitude of the soft compo- 
nent sees to be constant during the part of the decline when 
the spectral state is soft (Ebisawa 1991), indicating that the 
covering factor of the clumps cannot to be low, so that the 
disc effectively looks continuous. 

The highest luminosity spectra seem to need rather dif- 
ferent parameters of the Comptonization compared to the IS 
spectra. Perhaps, the disc structure changes when the disc 
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Figure 6. The surface density - mass accretion rate diagrams for 
an aPtot accretion disc around a 6 Mq black hole. The four thick 
curves are computed for R = 10, 30, 50 and 100 -R g . The parts of 
the curves with negative slope represent unstable solutions. The 
thin horizontal lines show estimated values of rh = M/Mrcdd f° r 
GS 1124-68, assuming d = 3.5 kpc. 



solution is located on the upper stable branch (so called 'slim 
discs': optically thick, with advective cooling; Abramow- 
icz et al. 1988). Alternatively, if the viscous heating scales 
with the gas pressure rather than the total (gas plus ra- 
diation) pressure then the instability is avoided (Stella & 
Rosner 1984) . While such a scaling might be expected for an 
MHD dynamo origin of the viscosity (Stella & Rosner 1984) , 
this has not yet been confirmed by numerical simulations. 
An interesting alternative viscosity prescription is implied 
by hydrodynamic experiments on rotating flows (Richard 
& Zahn 1999). For a Keplerian disc, the viscosity can be 
written as v = /3S1kT 2 , where r is the distance from the ro- 
tation axis and the coefficient (3 ~ 1/lZe, where Tie is the 
Reynolds number. This prescription gives disc solutions that 
are both thermally and viscously stable (Hure, Richard & 
Zahn 2001). 



6 CONCLUSIONS 

We model the spectra of Soft X-ray Transients in soft 
states using more physical models for the emission. We use 
reflection models as opposed to phenomenological broad 
line/smeared edge features, and Comptonization models 
which include the low energy cutoff at the seed photons as 
opposed to a continuous power law. We find that the strong 
soft component cannot be described by simple disc models 
(disc blackbody with colour temperature and/or relativis- 
tic corrections). There is additional broadening of the soft 
emission, indicating the presence of intermediate (between 
the cool main disc and the hot X-ray plasma) temperature 
material. This can be modelled either as an additional black- 
body, or as additional Comptonization. Both could be plau- 
sibly present. Magnetic reconnection above the disc could 
produce an X-ray heated spot underneath the flare, and/or 
the electrons in the flare itself may have a hybrid distri- 



bution, with a thermal tail at low energies which addition- 
ally Comptonizes the disc emission as well as a non-thermal 
power law to give the hard X-ray flux. Alternatively (ad- 
ditionally?) there could be a warm skin over the inner disc 
which Comptonizes the highest temperature emission from 
the disc. 

An ionized reflected component is always detected ir- 
respective of the detailed model used to describe the soft 
component. 
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